The Transit Light Curve Project. 
VIII. Six Occultations of the Exoplanet TrES-3 



Joshua N. Winn^, Matthew J. Holman^, Avi Shporer^, Jose Fernandez^, 
Tsevi Mazeh^, David W. Latham^, David Charbonneau^, Mark E. Everett^ 

ABSTRACT 

We present photometry of the exoplanet host star TrES-3 spanning six oc- 
cultations (secondary eclipses) of its giant planet. No flux decrements were de- 
tected, leading to 99%-confidence upper limits on the planet-to-star flux ratio of 
2.4 X 10"'*, 5.0 X 10~^, and 8.6 x 10~^ in the i, z, and R bands respectively. The 
corresponding upper limits on the planet's geometric albedo are 0.30, 0.62, and 
1.07. The upper limit in the i band rules out the presence of highly reflective 
clouds, and is only a factor of 2-3 above the predicted level of thermal radiation 
from the planet. 

Subject headings: planetary systems — stars: individual (TrES-3, GSC 03089- 
00929) — techniques: photometric 



1. Introduction 

The detection of reflected light from close-in giant planets is a difficult but worthwhile 
goal. Knowledge of the albedo is important because the stellar insolation is a critical factor 
in the atmospheric structure and the overall thermal balance of these planets (Guillot et 
al. 1996, Saumon et al. 1996, Seager & Sasselov 1998). Measurement of the planet's reflection 
spectrum, or at least the broad wavelength-dependence of the albedo, would provide clues 
about the dominant scattering mechanisms and constituents in the planetary atmosphere 
(see, e.g., Seager & Sasselov 1998, Marley et al. 1999, Sudarsky et al. 2000, Barman et 
al. 2001). The difficulty is that the reflected light is a very small fraction of the direct 
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starlight, producing a planet-to-star flux ratio at opposition of eA = px^Rp/a)"^, where Rp 
is the planetary radius, a is the orbital separation, and px is the wavelength-dependent 
geometric albedojl] Even for a planet as large as Jupiter, and an orbital separation as small 
as 0.05 AU, eA ~ 10""^ px- The reflected signal is also expected to vary over the orbital period, 
in a manner depending on the orbital inclination and the phase function of the planetary 
atmosphere. 

Early attempts to detect reflected light from close-in giant planets, by Charbonneau 
et al. (1999), Collier Cameron et al. (1999, 2002), and Leigh et al. (2003a,b), relied on 
optical spectroscopy. Those investigators sought the reflected copies of the stellar spectral 
lines by combining suitably Doppler-shifted spectra taken over a range of orbital phases. 
They did not detect the reflected signal and placed upper limit^ on the geometric albedo 
in the visual band of ~0. 1-0.5, subject to some ambiguity because of the unknown radius 
and orbital inclination of the planets they observed. Liu et al. (2008) recently revisited the 
spectroscopic technique and applied it to HD 209458b, for which the planetary radius and 
orbital inclination are known, but the available data provided only a weak constraint on the 
geometric albedo (0.8 ± 1.6 from 554—681 nm). 

Another detection method involves polarimetry. The Stokes parameters are expected 
to vary over the planetary orbital period because reflected light is preferentially polarized 
while direct starlight is unpolarized (see, e.g.. Hough et al. 2006). Berdyugina et al. (2007) 
reported a detection of a time-variable polarized signal from HD 189733. Taking the scat- 
tering radius to be the same as the optical radius measured through transit photometry, the 
implied geometric albedo is larger than 2/3, which is the geometric albedo of a perfectly dif- 
fusing sphere. While this is physically possible, and indeed some Solar system objects have 
geometric albedos exceeding unity due to strong backscattering, Berdyugina et al. (2007) 
preferred an interpretation in which the albedo is smaller and the scattering radius is larger 
than the optical radius measured during transits. 

For planets whose orbits are viewed close enough to edge-on that they undergo periodic 
occultations by the star, a powerful and conceptually simple method is available: when the 
planet is hidden by the star, the total light should diminish by the fraction eA. Edge-on 
systems are also advantageous because the planetary radius and orbital inclination can be 



^The geometric albedo is defined as the flux reflected by the planet when viewed at opposition (full 
phase), divided by the flux that would be reflected by a flat and perfectly diffusing surface with the same 
cross-sectional area as the planet. 

^Although Collier Cameron et al. (1999) reported a detection with greater than 95% confidence, additional 
data obtained by the same group did not confirm the detection; see Collier Cameron et al. (2002). 
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determined precisely from transit observationso If the measurement noise were limited only 
by photon statistics, the photometric technique could be employed with a smaller telescope 
than would be needed for either the spectroscopic or polarimetric techniques. However, 
ground-based photometry is generally afflicted by systematic errors that prevent one from 
achieving the required precision. To date, only spaceborne photometry has provided mean- 
ingful upper limits on exoplanetary albedos. In particular, Rowe et al. (2007) have used the 
Microvariahility and Oscillations of Stars satellite to set a 3a upper limit of 0.17 on the 
geometric albedo of HD 209458. 

Discoveries of transiting planets have abounded in the past few years, and there are 
now several systems known with larger values of (i?p/a)^. One of the most favorable systems 
is TrES-3, discovered by O'Donovan et al. (2007), for which [Rp/aY = 7.5 x 10~^. In this 
system, a planet with mass 2 Mjup and radius 1.3 -Rjup circles a G dwarf star with an orbital 
period of 31 hr. This paper relates our attempts to detect reflected light from the TrES-3 
planet using ground-based photometry with meter-class telescopes. In § 2, we describe our 
observations and data reduction procedures. In § 3, we present our results. In § 4, we place 
these results in the context of other observations and of theoretical expectations for the 
albedo and the thermal emission spectrum. 



2. Observations and Data Reduction 

We observed TrES-3 on 6 different nights when planetary occultations were expected, 
according to the ephemeris of O'Donovan et al. (2007). Assuming that the planetary orbit 
is circular, as expected for a planet with such a short orbital periocEl, the uncertainty in the 
predicted mid-occultation time was smaller than 3 minutes on each night. Observations on 
other nights were attempted, but only the six nights described below offered clear enough 
skies for high-precision photometry. 

On 2007 April 27, 2007 July 4, 2007 September 14, and 2008 March 12, we used the 
1.2 m telescope at the Fred L. Whipple Observatory (FLWO) on Mount Hopkins, Arizona. 
(Here and elsewhere, the quoted date refers to the UT date at the start of the night.) The 



■^It is possible, however, for a planet on an eccentric orbit to exhibit occultations without transits, or 
transits without occultations (see, e.g., Irwin et al. 2008). 

''^The timescale for tidal circularization is ~2 x 10^ yr, using Eqn. (9) of Rasio et al. (1996) with the 
system parameters of O'Donovan et al. (2007) and a tidal quality factor Q — 10^. This is shorter than the 
estimated age of the star, 5 x 10^ yr (Torres et al. 2008). The current radial velocity data are consistent 
with a circular orbit and give an upper limit on the eccentricity of 0.12 with 99% confidence. 
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detector was KeplerCam, a 4096^ CCD with a square field of view 23'1 on a side (Szentgyorgi 
et al. 2005). We binned the images 2x2, giving a scale of 0'/68 per binned pixel. We obtained 
repeated 60 s exposures over the course of 4-5 hr bracketing the predicted mid-occultation 
time. The dead time between exposures was lis. Autoguiding maintained the pointing to 
within 2-3 pixels throughout the observations. On the first night, 2007 April 27, wc observed 
through a Sloan i filter (RiO.7-0.85 /xm), as the star rose from an airmass of 1.4 to the 
meridian. The full-width at half-maximum (FWHM) of stellar images was about 4.5 pixels 
(3"). On the next two nights, 2007 July 4 and 2007 September 14, we observed through a 
Sloan z filter (fsO. 85-1.0 //m, with the red cutoff arising from the quantum efficiency of the 
detector). On both of those nights, the target star began near the meridian and set to an 
airmass near 2.0 throughout the observations, and the FWHM was approximately 3 pixels 
(2"). On the last night, 2008 March 12, we chose the Sloan i filter again and observed 
through an airmass ranging from 1.1 to 2.0. Conditions were nearly photometric, and the 
seeing was exceptionally good and stable, with a FWHM of 2.1 pixels (1'.'4). 

On UT 2007 August 26 and 2007 October 3, we used the 1 m telescope at Wise Ob- 
servatory in Israel. We used a Princeton Instruments 1340 x 1300 back-illuminated CCD, 
giving a field of view of 13^0 x 12'.6 and a pixel scale of 0'.'58. The readout and refresh time 
was 25 s, and we used exposure times ranging from 60 s to 160 s, depending on the airmass 
and seeing. The guider generally maintained the pointing to within a few pixels, but there 
was a glitch on each night that led to a 10-20 pixel offset. We observed through a Bessell 
R filter (RiO.55-0.7 //m). On August 26, we observed for 6 hr as the target star set from 
the meridian to an airmass of 2.8, and on October 3, we observed for 4 hr as the target set 
from airmass 1.1 to 2.5. In both cases, the data obtained through an airmass greater than 
2.0 proved to be much noisier than the rest of the data, and ultimately were not used in 
our analysis. On both nights the seeing steadily worsened with increasing airmass, from a 
FWHM of 2.5 pixels to 4.5 pixels (1'.'8 to 3^2). 

We used standard IRAF procedures for overscan correction, trimming, bias subtraction, 
and fiat-field division. We performed aperture photometry of TrES-3 and several nearby 
stars. The sum of the fluxes of the comparison stars was taken to be the comparison signal. 
The flux timeseries for TrES-3 was divided by the comparison signal, and then by a constant 
chosen to give a unit mean flux outside of the predicted occultation. We experimented with 
different choices for the aperture size, and different combinations of the comparison stars, 
aiming to minimize the standard deviation of the out-of-occultation portion of the TrES-3 
light curve. Not surprisingly, the best results were obtained when the aperture diameter was 
about twice as large as the FWHM of the stellar images, and when as many comparison stars 
as possible were used with mean fluxes ranging from 50-150% of the mean flux of TrES-3. 
For the FLWO data, 10-13 stars were used, and for the Wise data, 6-7 stars were used. The 
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Wise light curve had an abrupt jump at precisely the time when the guiding failed and the 
pointing changed; to correct for this, we set the mean out-of-occultation flux to be unity 
for each of the two segments of data. Finally, for all data sets, we performed an airmass 
correction. Using only the out-of-occultation data, wc fitted the light curve to an exponential 
function of airmass, divided by the best-fitting function, and renormalized the light curve to 
have unit mean fiux outside of the occultation. 

The final relative photometry is plotted in Figure 1, and are given in numerical form in 
Table 1. One of the transit light curves of O'Donovan et al. (2007) is also plotted in Figure 1, 
because the transit data were used in our model for the occultation data, as described in 
the next section. The duration of the occultation is expected to be nearly the same as the 
duration of the transit, since the orbit is presumed to be circular, as discussed further in § 4. 

Since no fiux decrement is obvious in any of the light curves, the results are upper limits 
on the planet-to-star fiux ratio in each bandpass. This makes it especially important to 
understand the characteristics of the noise. The out-of-occultation FLWO i light curves from 
2007 April 26 and 2008 March 12 have standard deviations of 9.1 x lO"'' and 9.0 x 10"^, 
respectively. The expected level of photon-counting (Poisson) noise is 8 x 10~^, and the 
expected level of scintillation noise is 3 x 10~^, according to Eqn. 1 of Young (1967). The 
quadrature sum of these terms is 8.5 x 10~^. The FLWO z light curves from 2007 July 14 
and 2007 September 14 have standard deviations of 1.8 x 10^"^ and 1.4 x 10^^, as compared 
to the expected level of about 1.2 x 10^'^ from Poisson and scintillation noise. Thus, most of 
the observed noise can be accounted for. The noise in the z-band light curves, in particular, 
is only 6-7% larger than the noise from these two expected sources. 

The out-of-occultation Wise data from 2007 August 26 and 2007 October 3 have stan- 
dard deviations of 2.2 x 10^'^ and 1.9 x 10^^. In both cases the expected level of photon- 
counting noise is about 10^'^ and the expected level of scintillation noise is about 3 x 10^^. 
Thus, the noise in the Wise data is about twice as large as one would expect from these two 
sources. Because of the jump that was observed in each light curve at the time of a large 
pixel shift, we suspect that at least some (and perhaps most) of the excess noise is due to 
an imperfect flat field. On other nights when the telescope lost tracking more frequently the 
noise level was indeed even higher, although those other nights were also characterized by 
poorer sky conditions, making it impossible to isolate the source of excess noise. 

Next we tested for any correlations with time or other external variables, and for sig- 
nificant non-Gaussianity. We found no significant correlation (a correlation coefficient <0.2) 
between the final relative fiux data and the airmass, the pixel position, or the measured 
shape parameters (FWHM, ellipticity, position angle) of stellar images. Of course, wc had 
already decorrelated against the airmass, and for the Wise data we had already removed 
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Fig. 1. — Photometry of TrES-3. The top panel shows the 2;-band transit photometry of 
O 'Donovan et ah (2007). The other 3 panels show our photometry spanning occultations 
(secondary eclipses). The dotted hnes show the expected start and end of the occultation. 



- 7- 



the effects of the two pointing glitclies, so tlie lack of correlations with those variables was 
expected. For each light curve, we also examined the histogram of flux values, the auto- 
correlation function, and the factor by which the noise level is reduced when the data are 
binned. Because the histograms are approximately Gaussian, the autocorrelations are gener- 
ally as small as would be expected for uncorrelated noise, and the noise falls approximately 
as 1/ V^, in what follows we treat the noise as Gaussian and uncorrelated. 

3. Determination of Upper Limits on the Albedo 

One might imagine measuring the planet-to-star flux ratio, e\, by finding the ratio of 
the mean flux during the occultation to the mean out-of-occultation flux. However, because 
the TrES-3 planet is on a nearly-grazing trajectory, the durations of the ingress and egress 
cannot be ignored. Furthermore, the uncertainty in the orbital inclination i introduces some 
uncertainty in the expected shape of the light curve, and the conversion from ex into the 
geometric albedo p\ requires knowledge of the ratio Rp/ a and its uncertainty. 

To take all of these factors into account, we simultaneously fitted a parameterized model 
to our i?, z, and 2;-band occultation photometry and the high-precision B and z-band transit 
photometry of O'Donovan et al. (2007). The model and the fitting method are similar 
to those we have employed in previous papers in this series (see, e.g., Winn et al. 2007, 
Holman et al. 2007). The model posits a circular orbit of a planet with radius Rp and a 
star with radius i?*, with orbital inclination orbital separation a, and period P. The flux 
ratio between the full disk of the star and fully-illuminated disk of the planet is e^, where 
k refers to the bandpass. When the planet is projected in front of the star, or vice versa, 
we use the analytic formulae of Mandel & Agol (2002) to compute the appropriate flux 
decrement. We assumed that the flux received from the planet does not vary appreciably 
over the 5 hr span of our occultation data, i.e., that the planet's phase function is constant 
within 30° of opposition. For the star, we assumed the limb-darkening law to be quadratic, 
with coefficientcl taken from Claret (2000, 2004). The model parameters were ej, e^, e^j, 
Rp/a, R^,/a, i, P, and Tc (a particular time of midtransit), along with the airmass correction 



5 Specifically, we used a = 0.2508, b = 0.3019 for the z-band data; and a = 0.7250, b = 0.0967 for the B- 
band data. These are based on interpolations for a star with Tcff = 5720 K, logy = 4.6, and solar metallicity. 
Southworth (2008) has found that the uncertainties in the transit parameters (the scaled radii of the star 
and planet, and the orbital inclination) are underestimated when the limb-darkening coefRcients are fixed 
at theoretical values. However, our results for the geometric albedo are not susceptible to this problem, as 
we confirmed by repeating the fit with a linear limb-darkening law and allowing the coefficients to be free 
parameters. 
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parameters for each light curve. 

For the transit hght curves, we adopted the same flux uncertainties as O'Donovan et 
al. (2007). For the occupation hght curves, based on the investigation of the photomet- 
ric errors described in the previous section, we took the standard deviation of the out-of- 
occultation flux to be the uncertainty in each flux measurement. To determine the best- 
fitting values of the parameters and their uncertainties, we used a Markov Chain Monte 
Carlo (MCMC) algorithm. This algorithm creates a random sequence of points in parame- 
ter space whose density is an approximation of the a posteriori probability distribution of 
the parameter values. The sequence is generated from an initial point by iterating a "jump 
function," which in our case was the addition of a Gaussian random number to a randomly- 
chosen parameter. If the new point has a lower than the previous point, the jump is 
executed; if not, the jump is executed with probability exp(— Ax^/2). We used an isotropic 
prior for the orbital inclination (uniform in cosi), and a Gaussian prior in the orbital period 
to enforce consistency with the value P = 1.30619 ± 0.00001 d, which was determined by 
O'Donovan et al. (2007) using the yearlong baseline between the TrES survey data and the 
high-precision light curves. We used uniform priors for all other parameters. The planet-to- 
star flux ratios were required to be positive. After creating several independent chains to 
verify that they all converged to the same region of parameter space, we created one long 
chain of 5 x 10^ points for our flnal results. 

In Table 2, we provide the results for the 99%-confidence upper limit on the planet- 
to-star flux ratio and the geometric albedo for each bandpass. The geometric albedo was 
computed as p\ = t\{a/ Rp)^. The upper limit was defined as the value for which the cumula- 
tive a posteriori probability distribution took the value 0.99. Figure 2 shows the a posteriori 
probability distribution for the geometric albedo based on the data from each bandpass, and 
for all of the data in combination. The results for the transit-based parameters Rp/ B^, R^^/a, 
i, and Tc were in agreement with the published values of O'Donovan et al. (2007). 

4. Discussion and Summary 

This study has resulted in empirical upper limits on the planet-to-star flux ratio, and 
the geometric albedo, in the wavelength range 0.6—1.0 //m. The most stringent upper limit 
was obtained in the Sloan i band, where Pi < 0.30 with 99% confldence. To put this result 
in perspective, we remind the reader that a perfectly-diffusing sphere has a geometric albedo 
of 2/3. The TrES-3 planet is less reflective than such a sphere, at least for observations in 
the i and z bands. The bandpass-averaged geometric albedos of Jupiter in the i?, and z 
bands are approximately 0.50, 0.45, and 0.25, respectively (Karkoschka 1994). Our results 
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Table 1. Photometry of TrES-3 



Telescope Filter 


Heliocentric Julian Date 


Relative flux 


FLWO i 


2454217.806989 


0.99905 


FLWO i 


2454217.807834 


1.00095 


FLWO i 


2454217.808667 


1.00006 


FLWO i 


2454217.809524 


1.00156 


FLWO i 


2454217.810369 


1.00144 


FLWO i 


2454217.811214 


0.99904 



Note. — The time stamps represent the Hehocentric Juhan 
Date at the time of mid-exposure. We intend for this Table to 
appear in entirety in the electronic version of the journal. An 
excerpt is shown here to illustrate its format. The data are also 
available from the authors upon request. 



Table 2. Results of Occultation Photometry of TrES-3 



Bandpass 


Upper Limit (99% confidence) 




Planet-to-star fiux ratio Geometric albedo 


i 


0.00024 0.30 


z 


0.00050 0.62 


R 


0.00086 1.07 



-10- 




Fig. 2. — Constraints on the geometric albedo of exoplanet TrES-3, averaged over various 
bandpasses. Each curve shows the a posteriori probabihty distribution determined with the 
MCMC algorithm. 

show that TrES-3 is less reflective than Jupiter in the i band, although they do not enable 
a meaningful comparison in the other two bands. Firm comparisons to other exoplanets are 
also difficult at this stage because almost all published reports are upper limits. Using the 
spectroscopic method, Charbonneau et al. (1999), Collier Cameron et al. (2002), and Leigh 
et al. (2003a,b) have placed upper limits of 0.1-0.4 on various systems, although in those 
cases the results are subject to extra uncertainty because the planetary radius and orbital 
inclinations are not known. Rowe et al. (2007) used several months of spaceborne photometry 
of HD 209458 to conclude Popt < 0.17 for that system with 3a- (99.73%) confidence, over an 
optical bandpass ranging from 0.4-0.7 /^m. At the same confidence level of 99.73%, our 
upper limit in the Sloan i band is pi < 0.35. As mentioned in § 1, Berdyugina et al. (2007) 
found a surprisingly large planet-to-star flux ratio for HD 189733, indeed large enough that 
those authors suggested the scattering radius is larger than the transit radius. Our results 
do not permit such an interpretation for TrES-3. 

Theoretical reflection spectra for strongly irradiated giant planets such as TrES-3 de- 
pend upon many factors and are especially sensitive to the presence or absence of clouds. 
Without highly reflective clouds, they are generally expected to have very low optical albe- 
dos. The underlying cause of the low albedo in cloud-free models is strong absorption in the 
pressure- broadened resonance doublet transitions of sodium and potassium (see, e.g., Seager 
et al. 2000, Marley et al. 1999, Sudarsky et al. 2000) or, for the hottest planets, the plentiful 
electronic transitions of gaseous TiO and VO (Hubeny et al. 2003, Sharp & Burrows 2007, 
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Fortney et al. 2007). Condensates may form reflective clouds that can drastically raise the 
albedo. The particular condensates that are available will depend on the temperature and 
pressure profile of the atmosphere. 

As a specific example, Sudarsky et al. (2000, 2003) predicted that gas giants with 
effective temperatures between about 900 K and 1500 K ("Class IV" in their terminology) 
have optical albedos ^ 0.2 in the bandpasses we employed, while for temperatures ^ 1500 K 
("Class V"), the albedo is as large as 0.6. The difference is that the hotter planets have a 
high layer of opaque silicate and iron clouds in their upper atmospheres. The case of TrES-3 
is borderline, with an effective temperature of 1643 f^^^{l — AY^^ K (O'Donovan et al. 2007), 
where A is the Bond albedo (the ratio of reflected to incident bolometric power) and / is 
a phenomenological factor accounting for atmospheric circulation (with / = 1 for isotropic 
reradiation and / = 2 for radiation from the dayside only). For TrES-3 it would be possible 
to obtain a self-consistent Class V solution with highly-reflective silicate and iron clouds, 
e.g., / = 1 and A ^ 0.4, giving Tes = 1500 K. Our upper limit on pi rules this out, and 
hence may be regarded as evidence against the presence of highly reflective clouds in the 
planet's upper atmosphere. 

For cloud-free atmospheres, the thermal radiation from the planet is expected to produce 
a larger signal than the reflected light (see, e.g., Lopez- Morales & Seager 2007, Fortney et 
al. 2007). One might wonder whether our limits on the planet-to-star flux ratio also constrain 
the thermal emission spectrum of the planet. The thermally emitted flux can be estimated 
from the effective temperature of the planet. For TrES-3, even for the hottest possible 
temperature of 2100 K (for / = 2 and A = 0) we estimate that the maximum planet- 
to-star flux ratio due solely to thermal emission is 2 x 10"'^ at a wavelength of 0.9 /xm. 
This rough estimate agrees within a factor of 2 with a more detailed atmospheric model of 
TrES-3 by Fortney et al. (2007). Those investigators classify TrES-3 pM planefl for 
which they predict a very low albedo (because conditions are too hot for any condensates) 
and an "anomalously hot" brightness temperature (because of a high-altitude temperature 
inversion). For our observing bandpass the calculated planet-to-star flux ratio is ~10~^. 
This is smaller than our most constraining upper limit of 2.4 x 10~^ in the i-band, and hence 
we conclude that the expected level of thermal emission is beneath our detection limit by a 
factor of 2-3. 



^The distinction between pM and pL planets, in the nomenclature proposed by Fortney et al. (2007), 
is analogous to the M/L transition in low-mass stars and brown dwarfs. The pM planets are hot enough 
for TiO and VO to exist in gaseous form, leading to strong absorption of the stellar flux at low pressure, 
a temperature inversion in the planet's upper atmosphere, and a brightness temperature in excess of the 
equilibrium temperature in the optical and infrared bands. 
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In summary, through high-precision photometry of an especially favorable target, we 
have placed meaningful upper limits on the albedo of a transiting exoplanet. This is the 
first time this has been done with ground-based occultation photometry. Our results seem 
to be limited by random noise, suggesting that the acquisition of more data will lead to 
increased sensitivity. For example, with 8 more light curves comparable in quality to the 
two i-band light curves presented here, the photon-limited la error in the albedo would 
be approximately 0.03. Even if the albedo proves to be very small, as expected in cloud- 
free atmospheric models, such a data set would enable a ~4ct detection of the expected 
level of thermal emission. Of course, it is possible that sources of systematic noise will 
become limiting factors, such as flat-fielding errors, or time-variable differential extinction 
beyond a gradual airmass dependence. Nevertheless, our results give reason to hope for 
an unambiguous detection of reflected light or thermal emission from exoplanets based on 
ground-based photometry using meter-class telescopes. 

We thank Jonathan Fortney, Willie Torres, and Alex Sozzetti for helpful conversations 
and correspondence; John Southworth, for his code for computing limb- darkening coeffi- 
cients; and an anonymous referee for a timely and detailed review. This research was sup- 
ported by Grant No. 2006234 from the United States-Israel Binational Science Foundation 
(BSF), Jerusalem, Israel. KeplerCam was developed with partial support from the Kepler 
mission under NASA Cooperative Agreement NCC2-1390 (PI: D. Latham). 
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